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Abstract 



We discuss the process e^e^ W~W~ mediated by heavy Majorana neutrino exchange in the 
t— and u channel. In our model the cross section for this reaction is a function of the masses 
(mjv) of the heavy Majorana neutrinos and mixing parameters [UeN) originating from mixing 
between the ordinary left-handed standard model neutrinos and additional singlet right-handed 
neutrino fields. Taking into account the standard model background and contraints from low 
energy measurements, we present discovery limits in the {rriN, U^n) plane. We also discuss how 
to measure in principle the CP violating phases, i.e., the relative phases between the mixing 
parameters. 
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1 Introduction 



The question why the masses of the observed neutrinos are much hghter than those of the 
charged leptons is a central unsolved problem in particle physics. An attractive scenarios to 
understand this is the following: Adding right-handed singlet neutrino fields to the standard 
model, the resulting mass spectrum can be such that there are 3 essentially zero mass Majorana 
neutrinos and additional heavy Majorana neutrinos. While in the generic case the masses of 
these heavy neutrinos are too large to detect observable consequences at present and future 
colliders, it not excluded however, that their masses are in the range of a few TeV and that their 
coupling strengh to the charged leptons are rather large. For the latter case, the phenomenology 
for e+e~ collisions where heavy Majorana neutrinos can be directly produced via t- channel W 
and s- channel Z° exchange, has been worked out in detail in IQ. The main conclusion is that 
the cross sections and designed luminosities are large enough to find heavy Majorana neutrinos 
essentially up to the kinematical limit, i.e., up to m^r ^/s, provided that the relevant mixing 
angles are near the present bounds inferred from low energy data. A similar conclusion also 
holds for single Majorana neutrino production in ej collisions. Another attractive reaction to 
be discussed in detail in this paper is e~e~ — > W~W~, which is mediated by t- and u- channel 
Majorana neutrino exchange (see Fig. 1). As the neutrinos do not have to be produced in 
this reaction, one is potentally sensitive to neutrino masses which exceed ^/s. Therefore, the 
aim of this paper is to work out the discovery region for such neutrinos in the parameter space 
(masses and mixing angles), taking into account standard model background reactions. 

In the second part of this introduction we present the model and list the present bounds on 
masses and mixing angles inferred from low energy precision measurements. In section II we 
discuss the signal reaction e~e~ — > W~W~ , while section III deals with the simulation of the 
standard model background. In section IV we present the discovery limits in the {mN,U^pf) 
plane. Section V finally deals with the question of how one can measure in principle the relative 
phases of the mixing angles; the presence of such phases is a necessary condition to observe CP 
violation in reactions like e~e~ —>■ W~W~ . 

1.1 The Model 

As in the standard model (SM) only left-handed neutrino fields are present, a Dirac mass (m^i) 
term for the neutrinos cannot be written down; also a Majorana mass term [rriM) involving only 
the left-handed neutrino fields (ul) cannot be constructed consistently within the SM because 
this would require Higgs- Triplets which are absent in the SM. Consequently, neutrinos are 
exactly massless in the SM. There are, however, many extensions of the SM (like SO{10), Eq, ...) 
with k extra neutrino fields {ufj) which are singlets under the SM gauge group Gsm- In & generic 
case symmetry breaking then induces the mass term 

L~ = -TTZmoyR - ]^rnM + h.c , (1) 

where m/? is a (3 x fc) Dirac mass matrix and rriM is a (fc x k) Majorana mass matrix. As an 
example, the grand unified group SO(IO) can break to the SM group through the chain 

50(10) Gsm x U{1)b-l ^ Gsm ^ SU{3), x U{l)em , (2) 
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where Aqut, v' and v = 174 GeV are the respective breaking scales in this chain. In this 
example, the Dirac- and Majorana mass matrices are of the form 

mc, = g-v , niM = h ■ v' , (3) 

where g and h are matrices of Yukawa couplings. It is conceivable that the residual B — L 
breaking scale v' is as low as a few TeV] in this case also the Majorana mass term is expected 
to be relatively light. 

To get the mass eigenstates, one has to diagonalize the neutrino mass matrix by means 
of a unitary (3 + x 3 + /c) matrix U . The resulting mass eigenstates are 3 light Majorana 
neutrinos Vi with masses my. and k heavy Majorana neutrinos Ni with masses m^r-. Expressing 
the weak eigenstates in the charged current by the corresponding mass eigenstates 0, we get 
the W — e — N {W — e — v) vertices relevant for our process e~e~ W~W~: 



^WeuWeN 



(4) 



It is instructive to look for a moment at the simpler 1 + 1 dimensional case with only one 
left-handed and one right-handed neutrino field. The masses and the mixing angle are 
given by the see-saw formulae 

ml, mn 

m^ ^ , mjv ~ mM , UeN ~ • (5) 

rriM mM 

In this case, the bounds on the light neutrinos imply that the mixing angle UeN has to be very 
small; consequently the cross section in e~e~ W~W~ is much too small to be observed. 
In the general 3 + k dimensional case, however, it is possible, that the 3 light neutrinos are 
exactly massless while the mixing angles UeN remain finite (i.e. nonzero)! This means that 
there is some chance to find heavy neutrinos with rather large cross sections through the process 
e-e- W-W-. 

In the following, we assume the light neutrinos to be massless (minimal mass generating 
case) and the mixing angles UeN and heavy masses mN to be essentially free parameters. 



1.2 Constraints on masses and mixing angles 

Before listing the constraints from low energy measurements, we would like to point out a 
model constraint. Assuming that there are no Higgs triplet fields, the Majorana mass term for 
the left-handed neutrinos is absent (as in eq. (|I|)). This implies the relation 

J2^.uI + J2^nU!n = o , (6) 

u N 

or in the minimal mass generating case {m^, = 0) 

J2mNU',N = ■ (7) 

N 

We now list the experimental contraints: 

■^We assume without loss of generality that the mass matrix of the charged leptons has been diagonalized 
prior to the neutrino mass matrix. 
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• Charged current universality implies that |f/eArp < 4 x 10 ^. The LEP-1 data lead to 
similar bounds 0. 

• Heavy Majorana neutrinos mediate the rare decay /i ^ 67. Non-observation of this 
decay leads to the strong bound \ J2NUfj,NU*j^\ < 2 x 10^^. This bound, however, does 
not necessarily lead to a more stringent bound on UeN than charged current universality. 

• The non-observation of neutrinoless double beta decay PPou implies a bound on both, 
the light and heavy neutrino sector. The best present evidence on the non-observation of 
the (3(3ou reaction is from the experimental limit on the process ^^Ge — See~e~, with 
Ti/2 > 5 X 10^^ years. For light neutrinos the bound reads 

< few eV . (8) 

In the minimal mass generating case (mj, = 0) this bound is satisfied automatically. For 
heavy Majorana neutrinos the bound is of the form 

X |hadr. matrix element ^ < . (9) 

Tl/2 

A recent paper [Q shows that in the older literature the hadronic matrix elements have 
been strongly overestimated. Taking the new estimates for the hadronic matrix elements, 
the bound becomes 

< 6 X 10"^ TeV-^ . (10) 
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2 Signal 

In the minimal mass generating case, where the light neutrinos are exactly massless by defini- 
tion, only the heavy Majorana neutrinos {Ni) contribute to the process e^e^ W~W~ . The 
corresponding Feynman diagrams are shown in Fig. 1. For ^/s ^ m\Y the cross section is to- 
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Figure 1: Feynman diagrams for the process e^e^ ^ W . 



tally dominated by the production of longitudinally polarized W~. To leading order in mw/y/s 
and in tree approximation the angular distribution da{e~j^e~j^ —>■ W~W~)/d cos 9 is given by the 
expression |^ 
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d cos 9 Stt 
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(11) 



with t = — (s/2) (1 — cos 6') and u = — (s/2) (1 + cos 6'). The general relation between mass and 
mixing parameters given in eq. @ for all neutrino flavors ensures unitarity and renormalizabil- 
ity. As a consequence of eq. (|^ [or eq. (|)], the cross section in eq. (|Tl]) vanishes in the limit 
■\/s — > oo because the leading constant term becomes proportional to J2n ^nU^n- Another 
feature which is easily seen from eqs. (|1T]) and is that the heavy Majorana neutrinos have 
to be degenerate in mass in order to have a non-zero cross section. 

Note, that in the limit where the masses of all heavy Majorana neutrinos are much larger 
than ^/s, formula (|lT]) could be further approximated to read 

2 

, mw <^ <^ rriN ■ (12) 

In the subsequent numerical analysis the lightest of the heavy neutrino masses is singled out 
for simplicity in the scenario 

mw ^ Vs < rriN^ <^ mN2,'mN3, ■ ■ ■ ■ (13) 
The results are then parametrized by the mass mi^{= mi^^) and the corresponding mixing angle 

UeN- 

In the reaction e^^e^^ — Vr~Vr~the invariant mass of the W~W~ pair equals a/s; this is in 
contrast to all the background reactions we discuss later, where the invariant mass of the W~W~ 
pair reaches ^/i only in very special kinematical configurations. Therefore, this invariant mass 
is a essential variable to get rid of the background. To fully exploit this, we concentrate on the 
hadronic decay modes of the W~W~ boson pair where all the decay products can be detected 
(no neutrinos as in the purely leptonic or semileptonic decay mode). The branching ratio of 
the W~W~ pair into the three possible channel is given roughly by 

BR{W~W^ : both W decay hadronically) 

BR{W~W~ : both W decay leptonically) 

BR(W~W~ : one decays hadr., the other lept.) 

Before discussing the background it is instructive to see that the number of signal events 
(using eqs. (jl^) and (|T^ ) can be rather large: For y/s = 500 GeV, m^r = 1 TeV and f/f^ = 
4 X 10~^ the production cross section for W'W" followed by hadronic decay of both W~ is 
a = 1.2 fb; this corresponds to 60 events assuming a luminosity of L = 50 fb~^. 

3 Background 

The most important backgrounds come from the reactions e~e~ W~W~vv (1), e~e~ — *■ 
W~ Z^pe~ (2), e~e~ — >• Z^Z^e~e~ (3), and e~e~ — >• W~W~^e~e~ (4), where the gauge bosons 
decay hadronically and the charged leptons escape along the beam pipe. If the jet charges could 
be reconstructed, then of course only reaction (1) would remain as a backgound. In this study 
we assume that this is not possible, hence deriving rather conservative discovery limits. While 
(1), (2) and (3) have been calculated by Cuypers et al. 0], the process (4) has been included 



da 



r<2 



d cos 6 Stt 



E 

N 



TP 



36 

81 
9 

81 
36 

81 



4 



by Barger et al. As mentioned earlier, the invariant mass of the two vector bosons is in 
general much smaller than y/s, because of the additional leptons in the final state (see also 
Fig. 3). However, there is a tail in the invariant mass distributions of the gauge boson pairs 
in the background reactions extending up to ^/s, because the additional leptons are essentially 
massless. Given the fact that the cross section of the background reaction (4) is more than two 
orders of magnitude larger than the signal (using optimistic parameters for mjy and Uen) one 
should try to further reduce the background by imposing additional suitable cuts which do not 



affect the signal significantly. To do such studies, we wrote a Montecarlo simulation of the 
signal and the background. Reaction (4) is clearly the most serious background source, because 
first, it has by far the largest cross section and second, this cross section is dominated by 
configurations where the final state electrons disappear in the beampipe. We therefore only 
take into account this background reaction. We calculated reaction (4) in the Weizsacker- 
Williams approximation for left-handed electron beams, which means that we have to work out 
in a first step the cross sections for the subprocess 77 —>■ W~^W~ —>■ jets for polarized on-shell 
photons. By calculating directly the matrix elements with 4 fermions in the final state and then 
doing the narrow width approximation of the W^-propagators we have correctly incorporated 
the polarization effects of the intermediate vector bosons. In a second step we convolute the 
77 cross sections with the appropriate Weizsacker- Williams photon distribution functions. The 
corresponding three Feynman diagrams are shown in Fig. 2. As a check, we compared our 





Figure 2: Feynman diagrams for the background process e^e^ — W^W e e . Only those 
three diagrams are shown which contribute in the Weizsacker- Williams approximation. 

Weizsacker- Williams results for e^e^ W~W^e~ e~ with the exact results in ref. |^ and 
found good agreement. The resulting invariant mass distribution of the final state hadrons is 
shown in Fig. 3. Needless to say, this reaction is only a background, if the final state electrons 



disappear in the beam direction. However, as seen in Fig. 4 of ref. [jTO[ in these proceedings, 
this background can only be reduced by a few percent when imposing a realistic angular cut 
on the electron directions. 

In order to reduce the background we require the hadronic mass mhad to be larger than a 
certain critical value rrfy[^\. For each ^/s we do the analysis for the two typical values m'f[^\ = 
0.90 X y/s and m^"^ = 0.80 x ^/s, respectively. Figs. 4 and 5 illustrate that the typical average 
transverse momentum p±{aver.) of the four jets, defined as p±{aver.) = (1/4) [|p_|_|(l) -|- . . . -|- 
|p_l|(4)] is higher for the signal than for the background (after imposing the m^"^-cut). 

One therefore can further enhance the signal/background ratio by imposing a lower cut on 
the variable p±{aver.). 
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Figure 3: Invariant hadronic mass distribution in the process e^e^^ ^ e e W — e e + 
jets in the Weizsacker-Wilhams approximation. 

4 Discovery Limits 

To summarize, imposing cuts on the invariant hadronic mass, rrihad, and the average transverse 
momentum of the jets, p±{aver.), the background to the signal process e~[e~[ — > W~W~ — > jets, 
which is characterized by the two parameters ttln and UeN, can be significantly reduced. As 
discussed in section I, there are constraints from low energy data. While neutrinoless double 
beta decay is the most stingent constraint (see eq. (|T0|)) for neutrino masses below 660 GeV, 
the bound from charged current universality {U^]y < 4 x 10~^) dominates for larger masses. 
As in our plots the mass is mostly larger than 660 GeV we did not plot the bound from PPou- 
Requiring a signal/backgound ratio > 1 and the number of signal events to be > 20, we worked 
out the discovery limits in the {m^, ?7|^)-plane. The results for = 500 GeV are shown in 
Fig. 6, assuming a luminosity L = 50 fb~^ and requiring the average momentum of the jets 
p±{aver.) to be > 25 GeV. In this figure the horizontal line represents the bound from charged 
current universality and/or LEP-1 data. Imposing the cut on rrihad at 450 GeV the allowed 
region in the parameter space is above the solid curve and below the horizontal line. If rrihad 
is cut at 400 GeV the allowed region lies between the dashed curve and the horizontal line. 
From the figure one sees, that it is very crucial at which value the cut on rrihad is imposed. 
Cutting at 400 GeV, the allowed parameter space is essential empty; if the cut is imposed at 
450 GeV, one is sensitive to a rather large parameter space: masses up to 2 TeV could be 
discovered. In an ideal experimental situation, without any fiuctuations in ^/s, one would of 
course make the rrihad-cnt basically at ^/s. In the real situation this energy is smeared, however, 
by beamstrahlung and initial state radiation. Of course one could simulate these effects in order 
to find out the optimal value at which the rrihad-cnt should be imposed. Also the cut on the 
average transverse momentum was imposed in a rather arbitrary, but reasonable way; also this 
cut could be optimized in order to become sensitive to a somewhat larger parameter space. 
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Figure 4: Signal (solid line) and background (dashed line) as a function of the average transverse 
momentum p±{aver.) of the jets. The hadronic invariant mass is required to be rrihad > 450 
GeV and the other parameters are chosen to be -^i = 500 GeV, rriN = 1 TeV, U^j^ = 4 x 10^^. 

The corresponding results for ^/s = 1 TeV and L = 100 fb"^ are shown in Fig. 7. It is 
impressive that one is sensitive to neutrino masses up to 12 TeV, assuming the mixing angle 
at its present upper bound. 



5 CP violating phase from the interference region 



In this section we wish to discuss an alternative scenario 0, |7 
neutrinos according to 

In this case an interference can arise at the level of the single helicity amplitude e^e^ W^W^ . 
To make this more explicit, we rewrite the cross section (0) using the parametrization 



with two relatively light heavy 
• ■ (14) 
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qi = {niN^) ; ga = ("^TVa) ; <?! < ?2 , a, & > (real) . 
The differential cross section which does not depend on the overall phase A, then reads 
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1,21 



(15) 



(16) 



(17) 




Figure 5: Same as Fig. 4, but for rrihad > 400 GeV. 



In eq. ( [TTD c = cos 5 is the cosine of the relative phase between m^Vi U^Ni '^a^2 UIn2^ while 
the quantities t, u, s = t + u and X = X{t,u), are understood to be measurable at any suitable 
finite set of values for the kinematic variables. As seen from eq. (|17|), the reduced cross section 
X depends on the 5 real parameters 

a , b , qi , q2 > and c = cos 6 ; — 1 < c < 1 (18) 

The process e^e]^ —>■ W'W~ , in the sense of a leading approximation, does not exhibit direct 
CP (or T) violation. A nonvanishing CP violation only results, if in addition to the CP odd 
phase 6 also a CP even phase a, generated by the electromagnetic interaction, is involved 
(e.g. through Coulomb exchange in the initial state). This is in analogy to direct CP violation 
in weak decays of hadrons: S plays the role of a CKM angle, while a is to be identified with 
a phase generated by the strong interaction. In order to find CP violation experimentally, the 
pair of CP associated reactions e~[e~[ W~W~ and e^e^ — >• 1^+1^+needs to be observed with 
high statistics, an unrealistic endeavour for the two reactions at hand. High statistics and 
precision are (would be) mandatory because it is the difference of respective differential cross 
sections which reveals the sought CP asymmetry. This difference, which would be strictly zero 
neglecting electromagnetic initial- or final state interaction, is much smaller (by a factor a) 
than the individual interference patterns. 

As the interference pattern is much more readily observable - due to its leading character 
- than any actually CP violating effect, we concentrate on working out a strategy to etablish 
(or refute!) the CP violating angle 6, which enters the differential cross section through its 
cosine. Assuming the discovery of the reaction e^e^ — > W~W~at some initial cm. energy, 
predominantly due to the exchange of a single (noninterfering) heavy Majorana neutrino, it is 
well conceivable that the sought interference can be established when going to higher energies. 

From eq. (p^T]) it becomes clear that all the 5 parameters a, b, c, gi, and q2 (within their 
respective ranges as shown in eq. (0)) can be extracted in principle from suitable data. 
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Figure 6: Discovery limits in the {mi\f,U^^)-pla.iae for ^/s = 500 GeV and luminosity L = 
50fb~^. The horizontal solid line at C/f^ = 4 x 10~^ represents the upper bound on C/f^ from 
charged current universality. Requiring pj^{aver.) > 25 GeV and rrihad > 450 GeV, the allowed 
region in the parameter space is above the solid curve and below the horizonal line. Using the 
looser cuts p±{aver.) > 25 GeV and rrihad > 400, the allowed region lies above the dashed curve 
and below the horizontal line. 

assuming this interference scenario is realized in nature with a sufficiently large cross section. 
To emphasize again, a value for c = cos 5 different from ±1 would estabhsh the existence of CP 
violation. 

As this issue might become relevant in future, we decided to elaborate a detailed stategy to 
extract the values of the parameters from data. However, as the details of this discussion are 
necessarily technical, we relegate it to the appendix. 

6 Summary and Conclusions 

In this paper we have investigated what can be learnt about the (heavy) neutrino sector from 
the reaction e^e^ W^W^aX future colliders. In a first scenario we have assumed that only 
one of the heavy neutrinos is light enough to lead to a detectable cross section. For this 
case we have worked out discovery hmits in the {rriN: U^j^)-plan.e. We found that at the cm. 
energy ^/s — 500 GeV neutrinos with masses up to 2 TeV can be discovered, assuming the 
relevant mixing angle to be at its present bound originating from low-energy experiments. At 
^/s = 1 TcV, the discovery region is much larger: Neutrinos with masses up to 12 TcV can 
be discovered. For the discovery of heavy neutrinos e^e^ colliders are better suited than e+e~ 
machines, because in the latter case the neutrinos have to be produced; consequently, one is 
only sensitive to neutrino masses smaller than ^/s. However, once a heavy neutrino is found, 
its properties (like decay channels etc.) could be better investigated in an e'^e~ envirnment 
(with sufficient energy). Therefore, as in many other aspects, e+e" and e~e~ colliders are 
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Figure 7: As Fig. 6, but for -y/i = 1 TeV and luminosity L = 100/6^^. The solid line 
corresponds to rrihad > 900 GeV and p±{aver.) > 50 GeV, while the dashed hne is for rrihad > 
800 GeV and pj_{aver.) > 50 GeV. The horizontal line again represents the upper bound on 
t/gjv froKi charged current universality. 



complement ary. 

Finally, we have investigated a scenario where two of the heavy Majorana neutrinos are 
relatively light. There, in principle the relative phase 6 between the mixing parameters can be 
measured. This phase is responsible for manifest CP violation in e^e^ —>■ W~W~: If 5 7^ 0, vr, 
direct CP violation is established. For this reason we have discussed a rather detailed strategy 
how to extract this phase from data. 



Acknowledgments: We thank F. Cuypers, T. Han and C.A. Heusch for helpful discussions 
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Appendix 

A Stategy to extract the CP violating phase 

As seen from eq. the reduced cross section X is a rational function in the variables t and 
M, i.e., of the form X = g/h, where g and h are symmetric polynomials in t and u. For the 
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following it is useful to write this somewhat more explicitly as 

X = X{t,u)= f^'^^^'^'l , (19) 

where g^y = g^^ and h^p = hpa are the coefficients of the symmetric polynomials g and /i, 
respectively. In eq. (|T9|) the obvious summations are tacitly understood. These coefficients 
themselves depend on the set of parameters n = [a, 6, c, gi, 52] we finally want to extract, more 
precisely 

hap{qi,q2) , gf,u { a,b,c ; qi,q2) . (20) 

For mathematical reason, as we will discuss later, it is convenient to extract in a first step the 
coefficients g^^ and hap without making use of their functional dependence on the parameters. 
Only after having determined these coefficients we use the functional dependence in order 
to extract the parameter set n. This has the most welcome advantage, that the analysis 
done in this way becomes sensitive to scenarios where one heavy Majorana neutrino interferes 
with some "other mechanism" (with a dependence of the same form with respect to t and u) 
which also contributes to the process e^e^ W~W~ . As this process violates lepton number 
conservation, this "other mechanism" necessarily corresponds to new physics different from 
heavy Majorana neutrino exchange. 

A.l Extraction of the coefficients g^iy and hap 

From the explicit form of the functions /, in eq. (p!7|), one easily finds that the indices (=powers 
of t and u) /i, u, a, (3 in eq. (|19[) range from to 4. In general such symmetric polynomials 
have 15 independent coefficients. More detailed inspection of the polynomial h shows that the 
coefficient /i44 = 1, while the remaining coefficients of h depend on the parameters qi and q2- 
For the polynomial g the coefficients (700 = 0, goi = (and therefore also 5^10 = 0) while the 
other 13 coefficients depend on the parameters qi, q2, a, b, and c. In total there are 14+13=27 
nontrivial coefficients. We now put eq. ([19|) into the form 

h^p[ru''X]-g^,t>^u'' = . (21) 



In eq. (|2l|) a linearly factorized form is achieved, separating the nonredundant subset of deter- 



minable quantities 



Xap = ^{ + t^u'^ ) X ; T^, = i ( t^u'' + t^'u^ ) (22) 



from the set of quantities to be determined (=unknowns) 

(23) 

Taking into account the constraints on the coefficients g^,^ and hap mentioned in the text, eq. 
(|2T|) takes the inhomogeneous linear form 

^0 + X^p Kp + T^, {-g^,) = (^ A ^ 0) 

a/3 7^44 ; Xq = X44 = ( )^ X (24) 

^iv ^ 00,10,01 
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The expression A 0) in eq. (p^) reflects the fact that there are systematic deviations 
from the hypothetical relation in eq. ( plf ) - theoretical as well as experimental - beyond the 
statistical errors, where the latter are unproblematic from the mathematical point of view. 

We now form the two vectors of unknowns zi and associated determinables Xq, -D^, restrict- 
ing ourselves to the nonredundant set as specified in eq. (|23) 



{ } — { ; —g^u} ', {De} — { Xaj3 ', T^y] 

{ £ = 1, 14 } ^ { < a < /3 < 4 ; a/3 ^ 44 } (25) 
{^= 15,. ..,27} ^{0</i<z/<4; /iz/^00,01} 

The association Dg^ <^ { X^fs ; T^^^} for a 7^ /5 or /i 7^ z/ is 

f 2X^0 for a < p ^ i = l,...,U , , 

Di={ ^ ' ' (26) 

\lT^y for /i< z/^£= 15,...,27 ^ ^ 

in order to take into account the symmetric summation relative to the indices a ^ j3 and ^ v- 
Then eq. (|^) is cast into conventional linear form 

^0 + EEi D,z, = {]{^A^ 0) 

{Xo,D,,A, } = {Xo,D,,A, }(t,M) ^^^^ 

In eq. (pTf ) { , -D^ } are measured or exactly known functions, while A is an unknown function 
of the kinematic variables. The latter is set whenever possible, but should be kept in the 
equation in order to avoid eventual inconsistencies, which typically arise when the signal (or the 
full details of the assumed interference pattern) are insignificant relative to the systematically 
misinterpreted background, systematic and/or statistical measurement errors. 

To the nonredundant £ set, consisting of the functions and unknowns { , 2;^ , Xq ; A } we 
associate - one by one - a given {fc} set {k = 1, 27) in the following way: we choose a set of 
weight functions over the region of experimentally accessible kinematic variables t, u 

{wk} ■■ Wk = Wk {t , u) ; k = 1, ...,27 

Xo ^ ( Xo ) = J dtdu Wk Xq 

r (28) 

Dn ^ D ki = dtdu Wk De 



A ^ A k = J dtdu Wk A 

The weighting functions w include delta functions singhng out 

- individual points or idealized bins (t,u)k- 

- line integrals, e.g. angular integrals with fixed cm. energy. 

There exists an infinite number of different {k} sets, each one defined through a given choice 
of 27 weighting functions {wk } ; {k = 1, ...27). 
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We concentrate on a given set to be interpreted as sample case. Eq. (p7D is thus transformed 
to an inhomogeneous linear system of equations 

Xo + Dz = 0) , (29) 

where we use matrix notation in order to suppress component indices. 

In the present context the standard notions of linear algebra receive specific interpretations. 
The regular case applies to a nonvanishing determinant of D. Because of the approximate 
nature of eq. (^9[), a vanishing and approximately vanishing determinant are equivalent in a 
sense to be precisely defined. The irregular case corresponds necessarily to the situation where 
the data does not allow to derive the hypothetical interference pattern, i.e. where the sought 
parameter c = cos 6 cannot be significantly determined. While also this case constitutes useful 
information, we concentrate now on the more interesting regular case. 

In order to generate at least one {k} set, at least 27 bins are needed. We choose one more 
in order to enable error analysis. This presupposes that at least of order 140 to 240 events fall 
into these 28 bins. As an example we may choose the 28 bins by taking 4 fixed cm. energies 
and 7 angular bins for each; the actual choice, however, is irrelevant for the following. One 
then can form maximally 28 {k} sets leaving out one of the 28 bins in turn. We number each 
{k} set with the index {n} ranging from 1 to 28. 

Thus 28 matrices D^^^ result, which we individually subject to predetermined criteria f\ 
for qualifying as a regular matrix. These criteria will translate into lower bounds for the 
determinants depending on D^^^ 



Det L'^"> 



> Ci{n}) > . (30) 



In general not all matrices D^"'^ will pass the test in eq. (pUD. Thus the last condition defining 
the regular case is to specify the fraction of {k} sets for which the test in eq. ( ^OD is successful. 
To fix ideas we set this fraction to 75 %. 

Having finally defined the regular case, we assume in the following that the data qualifies as 
regular. We retain only those {k} sets, which pass the test in eq. (pO|). We (re)number them 
with the label {n+j 

n+ = 1,...,N+ ; 21<N+<28 . (31) 

The range for A^+, the number of regular {k} sets, corresponds to the 75 % limit required 
above. We compare the calculational effort required to establish the regular case with a direct 

minimalization procedure to determine the five parameters at hand. Following the present 
method the extraction of 28 x 27 x 27 + 28 = 20'400 quantities Xo,n, -D^"^ was required. This 
corresponds to between 7 and 8 values for each of the five parameters ( 7.3^ ~ 20'400 ). 

Now we determine the unknowns z from eq. (|29|) setting the systematic deviation vector 
A = for each retained {k} set 

{z = -D-^Xo)^^^y ; n+ = 1,...,N+ . (32) 

We can assign weights to each regular {k} set, in inverse proportion to the estimated quadratic 
standard deviation 

Pn ; EPn. = 1 • (33) 

•^We do not discuss here details of these criteria, straightforward to define, which are to be used to ehminate 
the irregular case. 
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In eq. (^) and in the following the index + for = 1, A^4- is omitted. We then form the 
mean value and relative error correlation matrix 

Z = T,Pn Z^n} ; oil' =T,Pn (34) 

Zi ze' 



The estimate of the error matrix gj^, in eq. ( P^ is only sufficient to estimate the correlation 
matrix for — 1 out of the N = 27 components zg. Nevertheless we can perform a 
minimization - not to be confused with a statistical - on the systematic deviations A„ in eq. 
(p9|), which we have set to zero in the first approximation step displayed in eq. (|32|). To this 
end we cast the set of eqs. (^2|) into the form 

(35) 



D^^^ ( z - Z{n} ) = A{„| ; = - ) ' Xt^ 



In eq. (^) we have used the euclidean norm A|^| = J2k{'^{n})l to evaluate the systematic 
deviation of A|„} from 0. 

We determine the minimum and standard deviation of the function (z) defined in eq. 

m 

x'-xLn=iz-{z)f (D^) iz-{z)) 



xLn =EPn [D^"'^ ( {z) - Z{n} ) 

(Z) = {DY' EPnDi-y^D^-} Z{n} 
{D^) = EPn D^""^ ^ D^"^ 



(36) 



In eq. (^6]) the symbol T denotes transposition. 

All quantities except the free variable z are functions of the determinables Xq"^ , D^"^, 
depending in a highly nonlinear way on the matrix elements of D. The regular case ensures 
that the (symmetric) matrix (-D^) is regular. 

The procedure outlined here is equivalent to a direct fit of the parameters obtained from 
the hypothesis in the form of eq. (1171), using the quadratic deviation function 



X"={X - \ae^'/'h + be'^'/'f^y (37) 



with one notable difference : the relation in eq. (pq), quadratic in the free variable z, is 
functionally exact. In the latter, hypothesis and the determination of determinables are fully 
separated. The functional dependence of the quantities 

on the arguments Xq „ , D^'^^ and p„ is known and exact. The first two arguments, defined in 
eqs. (^) - (^), involve, beyond the cross section, known auxiliary functions used to formulate 
the hypothesis but as such are independent of it. Their determination is subject to systematic 
and statistical measurement errors only. These can be analyzed independently to any test 
of the hypopthesis. 
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Furthermore all errors, including the incomplete theoretical knowledge of signal as well as 
remaining background amplitudes, are fully accounted for by the unknown quantities A|„} 
defining the function according to eq. (|35|). 

o contours 

The regular case ensures a pronounced minimum of Ax2(z) at 2; = (2;). This is not 
surprising, since the conditions in eq. ( pOf ) defining the regular case qualify the hypothesis as 
viable (necessary but not sufficient conditions). 

We consider the (/ — a) contours (/ = 1, 2, ...) bounding the regions of z values 

^X\z) = {z - l,z)f m {z-{z)) < pxLn (39) 

The probabilistic syntax inherent to the notion of a contour is only true under the following 
circumstances 

a) all systematic deviations are zero and the hypothesis is true. 

b) the systematic deviations A{„} are themselves independent stochastic variables and the 

hypothesis is true. 

Case a) above applies in an approximative sense to the siuation where statistical errors 
dominate. This will most probably be the case, whenever the lepton violating signals discussed 
here will become observable for the first time. 

Nevertheless the value of / to be applied to the following analysis is a matter of initially 
judicious choice to be followed by subsequent checks. The method at hand allows several such 
checks with the same data. The first consists in comparing the value for (z) with z defined in 
eq. (H 

Irrelevant systematic deviations imply the relations involving / in eq. (^) 

7^1 ; 7\Ln = 0{z^ g'z\\ (D^) II) (41) 

In eq. (0) g'^ denotes the variance matrix estimated in eq. (^) and || (-D^) || an appropriate 
norm, e.g || {D'^) \ \ = tr (-D^). 

The tests outlined above can be extended, using the same data, choosing new bins, thus 
enlarging the number of regular {k} sets. It is understood that throughout 

((^)-Z|„})2, {Z-Z{^yf « {zf ,Z^ (42) 

holds. If no satisfactory factor / = 0(1) can be found, this means that the hypothesis cannot 
be significantly established, without seperate understanding of systematic deviations. 

We continue the discussion on the contrary assumptions, i.e. that a satisfactory (/ — a) 
contour is determined with / ~ 1 such that 

AX^ ( Ztrue ) < P xLn ^ ^true ~ (z) (43) 
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A. 2 Establishing the interference scenario 



All disqualifying criteria being positively met, still does not prove the hypothesis of the inter- 
ference scenario defined in eqs. (14) - (0), because so far we have only determined successfully 



the vector 2;, which is nothing but the coefficients g^^i, and h^p in eqs. (^) and ( P5D [see also 
eq. (^)]. The vector z depends on the set of parameters n = [qi,q2,a,b,c\, in which we are 
ultimately interested in. For each point in the parameter set n we calculate z = z{ir). The 
so determined vectors z sweep out a 5 dimensional subspace, denoted by S5, when varying the 
five parameters n in the allowed range specified in eq. 



The (/ — a) allowed region in the parameter space n is the intersection of 5*5 with the z 



region bounded by the corresponding contour according to eq. (|39|) 



{f-a [U]} = {HI Ax' {z [R]) < PxLn} (44) 

The final test of our hypothesis of an interference scenario of two heavy neutrino fiavors, requires 
the region of allowed parameters (/ — cr [ n ]) not to be empty or concentrated very near to the 
bounding (/ — a) contour. The best values for the five parameters [ n ] are then determined 
by the minimum of the function Ax^ { z [ n ] ) in eq. 



AxM M n,,,, ]) = Mm[u]Ax' (z [U]) < p xLn (45) 

Amusingly, an interesting situation also arises if this test fails, whereby the specific hypothesis 
(i.e., the interference of 2 Majorana neutrinos) is definitely disproved. However, on the level 
of polynomials of the same degree as those defining the hypothesis, an alternative interference 
pattern is established by the consistent form of the (/ — cr) contour defined in eq. (^Hj). This 
positively would prove the interference of the one fiavor signal, considered as established ac- 
cording to the criteria discussed in the previous sections, albeit with an amplitude of another 
nature than the one considered, also violating lepton number! Interference with a standard 
model background amplitude is by definition excluded. 

We continue the discussion for the alternative case where the criteria in eqs. ( ^4]) and (|45|) 
are satisfied. The full amplitude for two interfering Majorana neutrinos can be reconstructed 
for all values of the parameters 

qi = (mN^)' ; q2 = {rriN^f ; gi < ^2 , a, 6 > 

(46) 

c ; -1 < c < 1 

in the allowed region { / — o" [ H ] }. 

On this basis additional tests of the consistency of parametrization with the data can be 
performed. Furthermore the absolute values of the mixing parameters | UeN^ \ and | UeN2 I iii 
eq. (|T5|) can be determined. These are subject to the constraints unrelated to the data sample 
considered here p 

I UeN, ? + I UeN, ^ < 4 X 10"^ , (47) 

eventually narrowing down the acceptable parameter region. 

Finally we turn to the most subtle of the five parameters, c, the cosine of the relative phase 
between rriM^ U^Ni U^N2- "^^^ significant determination of c relies on a solid error 
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analysis with respect to the other 4 parameters, since the error matrix must show a high degree 
of correlation between the errors of c and those of the other parameters. In the vein of the 
present discussion we assume all further tests to reveal no inconsistencies. Then we foresee the 
result 

-1 <C_ < Cbest < C+ < 1 (48) 



according to the (/ — o) criteria defined in eqs. (^31) , ( ^if ) and (^Sf) . 
Three distinct possibilities arise, in ascending order of interest: 

1) c+ is not significantly different - within (2 — a) say - from 1, nor is c_ from -1. 

2) either c+ is not significantly different from 1, or c_ insignificantly from -1. 

3) c+ is significantly different from 1, and c_ from -1. 

Case 1) is unlikely, given the severe criteria in particular those regarding the regular case, 
discussed above, but cannot be entirely excluded. The conclusion is that the interference 
pattern can be established for some part of the allowed parameter region but the error margin 
is too large for an actual determination of c. 

Case 2) establishes the interference pattern but does not exclude the CP conserving relative 
phases 5 = 0, vr. 

Finally case 3), which serves as motivation of this analysis, not only establishes - within the 
errors - the interference pattern but proves, albeit indirectly, the CP violating character of this 
interference. 

We shall conclude this discussion on the hypothesis that this last case prevails. 
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